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A simple kinematic model for the trajectories of Listeria monocytogenes is generalized to a dy- 
namical system rich enough to exhibit the resonant Hopf bifurcation structure of excitable media 
and simple enough to be studied geometrically. It is shown how L. monocytogenes trajectories and 
meandering spiral waves are organized by the same type of attracting set. 



Introduction. Listeria monocytogenes is a widely dis- 
tributed pathogenic bacteria which occasionally causes 
serious illness in humans. L. monocytogenes evades the 
host's immune system by living inside its cells. Proteins 
located on the surface of the rod shaped bacteria catalyze 
the polymerization of the infected cells' actin molecules 
and this activity propels the bacteria through the cy- 
toplasm [![. The underlying mechanism of actin-based 
motility is a subject of great interest both because L. 
monocytogenes is a deadly pathogen and because actin 
filament assembly plays a role in many forms of cell 
movement [2J . A useful feature of actin-based motility in 
L. monocytogenes is the "comet tail" of actin filaments 
which are left behind as a cell is transported Q. The 
"comet tails" provide a record of bacterial trajectories 
in the cytoplasm. These trajectories, which vary con- 
siderably between individuals, can be complicated and 
orderly at the same time. In the present work it is shown 
how these trajectories can be explained with a dynami- 
cal system that has a well known type of attracting set. 
The presence of this attracting set can account for the re- 
newal of actin-based motility after the bacterium divides 
and the persistence of motility as the bacterium invades 
neighboring host cells. 

L. monocytogenes trajectories are the result of a com- 
plicated interaction of proteins in a host cell's cytoplasm 
or in a cytoplasmic extract [I[ . There has been extensive 
research into the molecular mechanism underlying actin- 
based motility but they do not account for the compli- 
cated forms of L. monocytogenes' actin comet tails 0- 
[l2j . In 13 1 Shenoy et al. present a simple and remark- 
ably effective kinematic model for the trajectories of L. 
monocytogenes in a thin layer of cytoplasmic extract. 

In the Shenoy et al. model the effect of actin poly- 
merization is assumed to produce a net force on the cell 
body which points slightly off center and which causes 
the bacteria to spin about its long axis as it travels in 
two dimensions. The bacterium's velocity is given by a 
vector whose direction varies sinusoidally with time and 
whose magnitude is fixed. Choosing units of measure so 
the speed is 1, letting s stand for arc length, and 9 stand 
for the velocity's direction the Shenoy et al. model, in 
a non-dimensionalized form, is d6/ds = Slcos(s) where 
O > represents the maximum deflection from forward 
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FIG. 1: Six curves in the (x, y)-plane determined by d6/ds = 
k + cos(s). The inset for each (x, y) curve points to its cor- 
responding parameter values (Q,k). Each (x,y) curve starts 
at the point (0, 0) (marked by an open circle) in the direction 
9 — 0. In two cases with k = the curves exhibit linear drift 
(which follows the dotted lines). The paths are qualitatively 
the same for k, — 1/20 but show an overall tendency to veer 
from a straight course (as indicated by the dotted curves). 
The region of parameter space shown here only contains a 
portion of the parameter values that can replicate L. mono- 
cytogenes trajectories. 



motion. 

Since d9/ds equals curvature the non-dimensionalized 
form of the Shenoy et al. model gives a one parame- 
ter family of intrinsic equations for planar curves (i.e. a 
two dimensional analog for the Frenet-Serret equations) 
which exhibit qualitative changes as the parameter is 
varied. Shenoy et al. show that for small O the curve 
is sinusoidal, for f2 « 2.5 it resembles a figure eight, and 
for larger values of f2 the curve tends to turn succes- 
sively clockwise and counter-clockwise around a sequence 
of points (fig. [T]). There are many qualitatively different 
types of curves for values of f2 from to 16 and Shenoy 
et al. show that L. monocytogenes display most if not all 
of these types. 

Although this kinematic model can reproduce many 
of the complicated L. monocytogenes trajectories it does 
not explain how the forces on the cell body arise or ad- 
dress the stability of the motion. The question of how 
actin-based motility arises is important because the bac- 
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teria have to divide in order to proliferate inside the host 
organism and after cell division the protein catalyzing 
actin polymerization is redistributed on the bacterial sur- 
face 14, 151. The issue of stability is important because 
one of the functions served by actin-based motility is to 
enable L. monocytogenes to create "pseudopodal projec- 
tions" from one host cell into another and thus allow the 
bacteria to avoid the host's immune system [l|. If the 
dynamics underlying actin-based motility was not stable 
then obstacles in the bacteria's path could disrupt the 
bacteria's entry into the neighboring host cell. 

In the next section the model of Shenoy et al. is ex- 
tended to allow the angular displacement of L. mono- 
cytogenes trajectories to accumulate. The third section 
explains how this extension leads to a dynamical system 
with the same type of attracting set as seen with mean- 
dering spiral waves. The fourth section shows how the 
existence of the attracting set accounts for actin-based 
motility developing and persisting in L. monocytogenes. 

2. Generalizing the kinematic model. It is well 
known in ballistics that it is difficult for a projectile to 
travel in a perfectly straight direction. For a self pro- 
pelled bacterium in a viscous medium even a small asym- 
metry in the cell body can cause it to eventually deviate 
from a straight course 16[. On the other hand in the 



Shenoy et al. model the angular displacement remains 
within fixed bounds along the entire length of the curve 
and the overall trajectory conforms to a straight line de- 
spite the small scale oscillations in its direction. 

To improve the accuracy of the kinematic model it is 
worth considering previous studies on actin-based motil- 
ity. Rutenberg and Grant 17( related the overall cur- 



vature of the trajectories to the number of randomly lo- 
cated actin filaments propelling the cell. They treated 
the torque produced by the filaments as a constant for 
relatively long periods of time which led to trajectories 
with constant curvature, dO/ds = k. For trajectories 
with it / starting in the direction 6q = the angular 
displacement grows in proportion to arc length. 

Evidence for the secular dependence of angular dis- 
placement on arc length in the Rutenberg and Grant 
model was subsequently presented in the extensive ex- 
perimental study on actin-based motility in L. mono- 
cytogenes [HI. The study found bacteria trajectories 
that were nearly circular and whose angular displace- 
ments and path lengths were nearly proportional to the 
elapsed time. Consequently the angular displacements 
were nearly proportional to arc length. 

Here we combine the approaches of Shenoy et al. and 
Rutenberg and Grant into a single model and present a 
summary of the types of paths this model displays. The 
Shenoy-Rutenberg model is, in non-dimensional form, 
dO/ds = k + ficos(s) where Q, k are constants. For 
small k the paths are qualitatively the same as for k = 
but they veer from a straight course (fig. [T]). Shenoy et 



al. found it useful to modify their model in a few cases 
by adding a low frequency term but this still did not al- 
low the angular displacement to accumulate. Adding a 
constant term as we do here does allow the angular dis- 
placement to accumulate and it helps shed light on the 
effectiveness of the Shenoy et al. model. 

The Shenoy-Rutenberg model is comparable in form to 
a model by Friedrich and Jiilicher for the chemotaxis of 
sperm cells [19]. Both models determine the curvature of 
the path followed by the cells using a constant curvature 
term and a second term but they differ in the form of the 
second term. For the Friedrich and Jiilicher model the 
second term is a function of the chemoattractant concen- 
tration and the internal signaling network. The paths 
produced by the Friedrich and Jiilicher model depend on 
the form of the concentration field. 

It is not very difficult to determine the paths produced 
by the Shenoy-Rutenberg model. Let (x(s),y(s)) denote 
the arc length parameterization of a path. The addition 
of k leaves the model in the form of an intrinsic equation 
for planar curves so, for fixed values of the parameters, 
the solutions are congruent and we can focus on the initial 
condition (x,y,9) = (0,0,0). This gives 
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Changing the sign of either f2 or k yields congruent (x, y) 
curves so we can assume f2, k > 0. While the integral 
cannot be evaluated in terms of elementary functions the 
curves are symmetrical and made up of congruent copies 
of an arc of length %. Since the (x, y) curve is invariant 
under reflection about the y-axis we can reflect the arc 
for < s < 7r to obtain the arc for —it < s < ir. 

For non-integral k let r = cot(irk)x(n) + y(n). It can 
be shown that 



x(s + 2n) 
y(s + 2tt) — r 
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From this it follows that the (x, y) curve can be obtained 
by iteratively rotating the arc for — tt < s < tt about 
the point (0, r) which is the center of symmetry for the 
figure. 

For non-integral rational k = p/q {jp, q coprime) and 
il 7^ the (x, y) curve is closed with q-iold rotational 
symmetry. It is the union of 2q congruent arcs of length 
it. For irrational k and f2 ^ the (x, y) curve is quasiperi- 
odic in the plane. It is the union of an infinite number 
congruent arcs with length tt. 

For integer values of k we can think of r as having gone 
to infinity. It can be shown that y(s + 2tt) = y(s). To 
express the value of x at multiples of tt we can use the 
integral representation for Bessel functions 

i r 

Jz(—Q) = — / cos(kc7' + fi sin(<j)) da (3) 



(this integral representation does not apply to non- 
integer values of k). From this it follows that the (x,y) 
curve can be obtained by iteratively translating the arc 
for — 7r < s < 7r horizontally by the distance 2ttJ^( — SI). 
When —SI is a zero of the (x, y) curve is closed with 
length 2ir. Otherwise the (x,y) curve is the union of an 
infinite sequence of congruent arcs with length n. This 
generalizes a similar result from [13| for the k = case. 

The points of maximal curvature on an (x, y) curve 
occur where s is an even multiple of 7r, the points of 
minimal curvature occur where s is an odd multiple of 
7r, and the curvature varies monotonically in between. In 
a neighborhood of (0,0) an arc of the (x,y) curve lies 
above the horizontal tangent at (0, 0) and for k > SI the 
curvature is positive everywhere. 

For 1 = k > CI the (x, y) curve has the form of a 
trochoid with its "petals" lying in a row (fig. [2j. For 
small SI and < k < 1 the (x, y) curve has the form of a 
hypotrochoid with its "petals" on the outside. For small 
SI and 1 < k < 2 the (x, y) curve has the form of an 
epitrochoid with its "petals" on the inside. 

Flower like curves such as these are traced out by the 
tips of spiral waves propagating through excitable me- 
dia. Spiral waves occur in diverse systems with very dif- 
ferent underlying mechanisms. This includes aggregating 
myxobacteria which form macroscopic waves as cells glide 
across a two dimensional surface [22j . A spiral wave often 
propagates as though it were a rigid body rotating about 
a quiescent core. Away from the core the shape of the 
wave front converges to an Archimedean spiral U 24$ . 
However under appropriate circumstances the inner tip 
undergoes a secondary oscillation as the wave rotates and 
thereby traces a hypo/epi/trochoid like curve. Spiral tip 
meander has been observed in many systems such as the 
BZ chemical reaction [2^|, heart tissue [2f|, and aggre- 
gating cells of Dictyostelium discoideum (cellular slime 
molds) [13. 

3. Resonant Hopf bifurcations. An important 
step toward understanding why spiral tip meander oc- 
curs in systems with such different underlying mech- 
anisms was made by Barkley 2H-3 



who recognized, 
through numerical and mathematical analysis, the im- 
portant role played by the group of orientation preserv- 
ing congruences of the Euclidean plane, SE(2), and that 
this role can be exemplified by reducing the dynamics 
to five dimensions. The mathematics of Barkley's break- 



through has been further elaborated and generalized [31 



35] . Barkley's approach can be nicely illustrated with the 
Shenoy-Rutenberg model since it already has the form 
of an intrinsic equation for planar curves. To do this we 
couple the Shenoy-Rutenberg model to a two dimensional 
system from [36[ . We write the Cartesian coordinates for 
this subsystem as (X,Y). The full differential equation 
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FIG. 2: A version of a Zykov-Winfree flower garden [37l [3S| 
whose isogonal contours have been combed straight. The isog- 
onal contours of d6/ds = k + S7cos(s) for k = 2/3, 4/5, 1, 6/5 
are shown. Each [x, y) curve in the insets starts at the point 
(0, 0) (marked by an open circle) in the direction 6 — 0. For 
k — 1 the (x, y) curves exhibit linear drift. For k below 1 the 
(x,y) curves have hypotrochoid like shapes and for k above 
1 the (x,y) curves have epitrochoid like shapes. So long as 
k > SI (above the diagonal line) the (x,y) curves do not have 
inflection points. 
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where the parameter /i corresponds to Barkley's normal- 
ized bifurcation parameter [29(. For /i < the origin of 
the (X, Y) subsystem is an attracting fixed point. At 
fi = the Hopf bifurcation occurs and for /i > there 
is an attracting circular limit cycle centered at the origin 
with radius yfji. 

To concisely express the initial conditions we set STi = 
for \i < and SI = J~\i for [i > 0. For the initial 
condition (x,y,0,X,Y) = (0,0, 0,S7, 0) the solution to 
the (X,Y) subsystem is (X(s) 7 Y(s)) = Sl(cos(s), sin(s)) 
which gives 6' = k + SI cos(s) which in turn recovers 
eq. (fTJ) for the (x, y) subsystem. 

A purely rotating spiral wave appears motionless in 
a frame rotating with it. The transition to meandering 
corresponds to the Hopf bifurcation. After the bifurca- 
tion the spiral tip appears in the rotating frame to trace 



4 



a circularly shaped path although far from the core the 
wave continues to appear motionless. 

By converting eq. (|4]) to a rotating coordinate sys- 
tem (0, 0, 0, 0, 0) becomes a fixed point with spectrum 
{±ik, 0, /i±i}. The eigenvalues ±ik arise from the trans- 
lational symmetry of the plane and arises from the ro- 
tational symmetry of the plane. At the Hopf bifurcation 
all five eigenvalues lie on the imaginary axis. 

Barkley showed that the type of curve traced by a spi- 
ral tip in the stationary frame depends on where the Hopf 
eigenvalues cross the imaginary axis in relation to the 
translational eigenvalues. When the translational eigen- 
values are between the Hopf eigenvalues the spiral tip 
will follow a hypotrochoid like curve (0 < k < 1 in 
eq. j4|). When the translational eigenvalues are outside 
of the Hopf eigenvalues (but not more than twice the 
Hopf eigenvalues) the spiral tip will follow an epitrochoid 
like curve (1 < k < 2 in eq. Q). When the translational 
and Hopf eigenvalues coincide the spiral tip exhibits lin- 
ear drift (k = 1 in eq. dU). 

In terms of L. monocytogenes we can interpret {X, Y) 
as the projection of the cell's translational velocity to 
a plane orthogonal to the cell body's long axis and we 
can interpret the oscillation of (X, Y) as the effect of 
the cell's spin on its propulsion system. The long axis 
and the X component are parallel to the surface being 
traversed while the Y component points in the orthogonal 
direction. For a cell constrained in two dimensions the Y 
component does not contribute to the motion. For Q = 
the cell appears motionless in a frame rotating with it. 
For small fi > the cell appears to follow a circularly 
shaped path in the rotating frame. 

The detailed mechanisms behind spiral meander and 
L. monocytogenes motility arc different but the paths 
they follow are both part of a larger two parameter fam- 
ily of curves. The paths followed by spiral wave tips 
are organized around a first order resonant Hopf bifur- 
cation for which the translational and Hopf eigenvalues 
coincide (k = 1 in eq. (jJJ). The paths followed by L. 
monocytogenes are organized around a zero order res- 
onant Hopf bifurcation for which the translational and 
rotational eigenvalues coincide (k = in eq. Q). 

4. The dynamics of L. monocytogenes motility. 

Spiral waves appear in excitable media when, in the state 
space for the medium, the state is sufficiently close to 
the appropriate attracting set. Each state in the attract- 
ing set corresponds to a well formed spiral wave in the 
medium. In many cases there is a characteristic wave- 
length to the limiting form of the spiral wave [3{|. In 
such cases any two spiral waves in a planar homogenous 
isotropic medium will be congruent. For media which 
support non-meandering spiral waves the attracting set 
is essentially a copy of the symmetry group SE(2). This 
is the type of attracting set the solutions to eq. ((H) have 
when fi < 0. The orbits of the dynamical system in- 



side the attracting set are simple closed curves and thus 
bounded. These orbits correspond to spiral waves un- 
dergoing a pure rotation or L. monocytogenes following 
a circular trajectory. 

Aside from numerical simulations it is difficult to pre- 
pare excitable media so that the initial state of the sys- 
tem is within the attracting set, i.e. so that the medium 
begins with a well formed spiral wave which then under- 
goes a pure rotation. A purely rotating spiral wave only 
appears after a transient period. One way for a purely 
rotating spiral wave to appear is by disrupting a circular 
or linear wave front with an obstacle in the medium. The 
broken end of the wave front will then curl up and over 
time the shape of the wave will develop into a well formed 
spiral which propagates in a purely rotational manner. 
The disruption of a wave front by an obstacle brings the 
state of the system sufficiently close to the attracting set 
that it converges towards it. 

For homogenous isotropic excitable media which sup- 
port meandering spiral waves the attracting set has an- 
other dimension. This is the type of attracting set the 
solutions to eq. ((H) have when /i > 0. Each point in 
the attracting set gives the position and orientation of 
the spiral wave as well as its phase within the period of 
meander. The orbits of the dynamical system inside the 
attracting set are bounded unless there is a resonance be- 
tween the rotation of the spiral wave and the oscillation 
of the tip in which case they are unbounded. Meander- 
ing spiral waves appear after a transient period once the 
state of the system has been brought sufficiently close to 
the attracting set. For systems at or near resonance the 
core of the spiral wave will be transported across large 
distances. 

Biological systems repeat many of the same develop- 
mental strategies in various contexts to form functional 
patterns. The presence of low dimensional attracting 
sets in complicated dynamical systems can provide sta- 
bility to developmental processes which are exposed to 
the environment. For instance there are prokaryotes 
{e.g. myxobacteria) and eukaryotes {e.g. D. discoideum) 
which use spiral wave dynamics to get individual cells 
dispersed over a wide area to aggregate together and de- 
velop multicellular reproductive organs. The underlying 
mechanisms by which myxobacteria and D. discoideum 
move and communicate are quite different but the pres- 
ence of an attracting set for spiral wave dynamics can 
allow the aggregation process to proceed despite the va- 
garies of their environments. 

There has been a long running and continuing effort to 
determine the underlying mechanism responsible for the 
formation of actin-based motility 0- 12 , 4(| . The Shenoy 
et at model is effective at duplicating L. monocytogenes 
trajectories but it is not directly based on a physicochem- 
ical mechanism. Their model proceeds from general con- 
siderations about how the forces produced by actin poly- 
merization must act on the cell. In order for the cell to 
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change direction as it moves there must be some asymme- 
try in the distribution offerees exerted on the cell surface. 
By treating the net propulsive force as a constant paral- 
lel to the long axis of the cell and whose exertion point 
rotates at a constant distance about the long axis the 
magnitude of the component of the net torque orthogo- 
nal to the plane of motion varies in a precisely sinusoidal 
fashion. In this way the cell body oscillates about its 
center of mass much like an ideal torsional spring. With 
the propulsive force always parallel to the long axis the 
cell moves in trajectories that alternately wind clockwise 
and counter-clockwise. 

However this clockwork like mechanism does not sim- 
ply appear fully formed. For L. monocytogenes engaged 
in actin-based motility the concentration of ActA (the 
catalyst for actin polymerization on L. monocytogenes) 
along the cell wall increases from the apical pole to the 
basal pole. L. monocytogenes cells reproduce by dividing 
along a septum midway between the poles. After divi- 
sion each daughter cell forms a new apical pole at the 
septation region and ActA concentration is redistributed 
along the cell wall 14, 15j |. 

When a L. monocytogenes cell hrst begins moving it en- 
gages in a "hopping" type motion. The density of actin 
builds up behind the basal pole until there is a sudden 
acceleration of the cell. The actin tail then becomes rar- 
efied, the cell subsequently decelerates nearly to rest, and 
the cycle repeats. After several cycles the cell eventually 
settles down to a relatively constant speed @, EH ■ 

One function served by the actin-based motility of L. 
monocytogenes is the transport of bacteria from one host 
cell to another without the bacteria having to leave the 
confines of the host cells. This is accomplished when a 
bacterium presses against the host's plasma membrane to 
create a "pseudopodal projection" with the bacterium in- 
side. The bacterium enters a neighboring host cell when 
the pseudopodal projection is phagocytosed by the neigh- 
boring cell 

During the transient period L. monocytogenes move- 
ment is sensitive to obstacles in its environment but it 
becomes more robust when it reaches a steady speed 
The presence of a low dimensional attracting set 
for actin-based motility can provide stability in the de- 
velopment of L. monocytogenes infectiousness. When the 
state of the system is in its transient phase, away from the 
attracting set, obstacles in the environment can have the 
effect of perturbing the system from one orbit to another 
with a very different course. On the other hand when 
the state of the system is close to the attracting set it 
can quickly return to the attracting set after an obstacle 
causes a perturbation. This accounts for how L. mono- 
cytogenes motility can persist as it forms a pseudopodal 
projection. 

Moreover the attracting set for actin-based motility by 
L. monocytogenes confined to two dimensions appears to 
be of the same type as for spiral tip meander in two di- 



mensions. In both cases the points of the attractor cor- 
respond to the position, the orientation, and the phase 
in the secondary oscillation of the object that is moving. 
For the purposes of transport it is useful for the system 
to be near a resonance. The Shenoy et al. model corre- 
sponds to a zero order resonance which helps account for 
its effectiveness although the model needs to be placed in 
a larger mathematical context to see this and to account 
for how actin-based motility in L. monocytogenes arises 
and persists in the presence of obstacles. 

Conclusion. Spiral waves involve the cooperative be- 
havior of multiple agents while a bacterium is generally 
thought of as a single agent. However this apparent asso- 
ciation of the first order resonant Hopf bifurcation to the 
organized motion of multiple agents and the zero order 
resonant Hopf bifurcation to the motion of a single agent 
is not a general principal. The movement of an individual 
bacterium can be regarded as the action of a single agent 
but actin-based motility involves the polymerization of 
actin so it can also be regarded as an organized activity 
involving many chemical agents. Inert beads coated with 
proteins which catalyze actin polymerization can also en- 
gage in actin-based motility 4l|, |42 1 . 

Actin is an important constituent of the cytoskeletons 
of eukaryotic cells and it forms locomotory structures 
such as filopodia in which actin organizes into bundles 
and lamellipodia in which actin organizes into meshworks 



49(. Numerical simulations indicate that actin can 
also form spiral waves [13] and there is evidence for actin 
for ming spiral waves inside of D. discoideum pseudopodia 
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The actin rich cytoplasm of eukaryotic cells can be 
seen as a type of excitable media with a propensity to 
locomotory behavior which several pathogens, in addition 
to L. monocytogenes, have evolved to take advantage of 
in order to transport themselves. Actin-based motility 
occurs with the bacteria Shigella flexneri [53| , species of 
Rickettsia bacteria 54 1 , and vaccinia virus particles (HBJ . 
The detailed mechanism of actin-based motility varies 
between these pathogens but the y ea ch involve catalyzing 
the polymerization of actin, [56l-l58j . 

When L. monocytogenes bacteria are constrained to 
move in two dimensions their trajectories can form a com- 
plicated series of coils. As with the trajectories of mean- 
dering spiral waves, the trajectories of L. monocytogenes 
are composed of essentially congruent copies of a finite 
curve repeatedly joined end to end. The form of the re- 
peating unit appears after a transient period and it varies 
between different occurrences of meandering spiral waves 
and between different individual L. monocytogenes. Yet 
complicated curves characteristic of both types of phe- 
nomena can be replicated using a two parameter family 
of dynamical systems with the same type of attracting set 
as eq. ([!]) . Taken together these lines of evidence support 
the idea that actin-based motility in L. monocytogenes is 
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organized by the same type of low dimensional attracting 
set that organizes spiral tip meander. 
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